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Abstract 

We consider the notion of thermal equilibrium for an individual closed macro¬ 
scopic quantum system in a pure state, i.e., described by a wave function. The 
macroscopic properties in thermal equilibrium of such a system, determined by its 
wave function, must be the same as those obtained from thermodynamics, e.g., 
spatial uniformity of temperature and chemical potential. When this is true we 
say that the system is in macroscopic thermal equilibrium (MATE). Such a sys¬ 
tem may however not be in microscopic thermal equilibrium (MITE). The latter 
requires that the reduced density matrices of small subsystems be close to those 
obtained from the microcanonical, equivalently the canonical, ensemble for the 
whole system. The distinction between MITE and MATE is particularly relevant 
for systems with many-body localization (MBL) for which the energy eigenfuc- 
tions fail to be in MITE while necessarily most of them, but not all, are in MATE. 

We note however that for generic macroscopic systems, including those with MBL, 
most wave functions in an energy shell are in both MATE and MITE. Eor a classi¬ 
cal macroscopic system, MATE holds for most phase points on the energy surface, 
but MITE fails to hold for any phase point. 
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1 Introduction 


Thermal behavior of closed macroscopic systems in pure states has been widely studied in 
recent years; see e.g., [H [30l UHl SSI |3Tl |32l |35l |23l HH |33l [36l |3ll [8] , after some pioneering 
work even earlier [38l 061 |39l |6l 001 01]. In particular, the importance of the eigenstate 
thermalization hypothesis (ETH) [2100] has become widely appreciated, see e.g., [3211321 
[2310210I112210111H]. It asserts, in one version, that all energy eigenstates (in a suitable 
energy shell) are thermal, i.e., assigning probability distributions to observables that are 
characteristic of thermal equilibrium. The ETH holds for many (but not all) macroscopic 
quantum systems. When it holds then the system, starting out of equilibrium, will 
thermalize (at least in the time average; see Section 0] below). Thus if a system does 
not thermalize it must have energy eigenfunctions that fail to be thermal. 

An important case of this is that of many-body localization (MBL) [H 01 02], for 
which the Hamiltonian has (at least some) eigenfunctions that are in some way localized, 
so that, for any wave function, the component from these eigenfunctions will not spread 
but stay localized forever. For systems with MBL it has been argued that most (if not 
all) energy eigenfunctions (in suitable energy intervals) fail to be thermal, and there are 
in fact models for which this can be analytically [12], numerically [02], or perturbatively 
[21 HI [37] seen to be the case. At the same time it has been argued [15] that most energy 
eigenstates must rather generally be thermal, in particular even for systems with MBL. 

This situation is thus rather puzzling: How can the fact that most energy eigen¬ 
functions for (some) MBL systems are not thermal be reconciled with the results in [T5] 
showing that in wide generality most of them must be. The answer to this question, 
we point out, lies in the fact that there are basically two notions of thermal equilibrium 
for pure states: a macroscopic notion of thermal equilibrium that we call MATE, and 
a more refined microscopic one that we call MITE. While most, but definitely not all, 
energy eigenstates of a system with MBL are in MATE (see Sec. 0]), none, or nearly 
none, are in MITE [21 El ESI 0S]. Nonetheless, most pure states in the energy shell are 
in both MATE and MITE, even for systems with MBL (see Sec. (21 and ED- It is with 
these two notions of thermal equilibrium and their consequences that we are concerned 
here. 

To be more precise, consider a finite, macroscopic, closed quantum system with 
Hilbert space Jif. Let be a micro-canonical energy shell, i.e., the subspace of ^ 
spanned by the energy eigenstates with eigenvalue in an energy interval that is small 
on the macroscopic scale but contains many eigenvalues. The micro-canonical density 
matrix is defined by = (dim=^^c)~^-i7nc with Pmc the projection to As 

usual, pure states in M'-iac are superpositions of energy eigenstates in Both MITE 

and MATE can be expressed as subsets of the unit sphere in 

S(^mc) = {V' e JCe : m = 1} . (1) 

We will often simply say “micro” for “microscopic” and “macro” for “macroscopic”; 
“most” means “all but a few” or “all except a set of small measure” (i.e., “an over¬ 
whelming majority of”); measures are taken to be normalized; the small measure, in 
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fact, tends to 0 in the thermodynamic limit. We dehne MATE and MITE in Sections [2] 
and [31 respectively, and discuss the properties and differences of the two notions in Sec¬ 
tion jH We also note there that these notions can also be applied to mixed states. We 
provide a deeper and more detailed discussion elsewhere [13] . 

2 MATE 


The dehnition of MATE is based on macro observables Mi,..., Mk- To be specihc, 
macro observables can reasonably be based on a partition of the system’s available 
volume A C into cells A* that are small on the macro scale but still large enough 
to each contain a large number of degrees of freedom. Examples of natural choices of 
M’s are, for each cell, the number of particles of each type, the total energy, the total 
momentum, and/or the total magnetization. These are the variables usually considered 
in the thermodynamic or hydrodynamic description of a fluid or magnet. 

Following von Neumann |l6], we take the Mj to commute with each other and to be 
such that the gaps between the eigenvalues are of the order of the macroscopic resolu¬ 
tion (so that the eigenvalues are highly degenerate). This can be achieved by suitably 
“rounding off” and coarse-graining the operators representing the macro observables 
[36l[I7l|27|. Taking to be an eigenspace of a “macro energy” operator, and thus 
to commute with the other macro observables, all Mj can be regarded as operators on 
^rac- Their joint spectral decomposition dehnes an orthogonal decomposition 

Mfmc = , 

V 

and the subspaces M'y (“macro spaces”), the joint eigenspaces of the macro observables, 
correspond to the different macro states and have very high dimension. The decom¬ 
position is in some ways analogous to a partition, in classical mechanics, of an energy 
shell Tmc in phase space into disjoint subsets Vy corresponding to different macro states 
[311H11I1E2]. It is generally the case llEQill], both in classical and quantum mechan¬ 
ics, that one of the macro states, corresponding to thermal equilibrium, is dominant, 
i.e., that one of the r,^’s, denoted Teq, has most of the phase space volume of Tmc, and 
that one of the Mi, denoted ^q, has most of the dimensions of i-e., 

dim J^q 
dim 

with e <C 10 Realistic values of e, say for a liter of air under atmospheric pressure, 
are smaller than 10“^° [13]; more generally, e is exponentially small in the number of 

^An exception to the existence of a dominant macro space is provided by first-order phase transitions, 
such as in the ferromagnetic Ising model in a vanishing external magnetic held, where Mi has the 
appropriate majority of spins up and Mi' has the appropriate majority of spins down, each having 
nearly 50% of the dimension of Mine for a suitable energy interval. 
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degrees of freedom per cell. We assume here that our system and our choice of macro 
observables are such that ([3]) holds for suitably small e. 

The system is said to be in MATE whenever its wave function -0 lies within a 6- 
neighborhood of Jifgq with suitably small 5 > 0, i.e., in the set 

MATE = |0 e §(JCc) : (01 Aq|0) > 1 - h} (4) 

with Peq the projection to J^eq- Thus for a state 0 that is in MATE, the probability 
is close to one that all macro observables take on their thermal equilibrium values. A 
concept of thermal equilibrium along these lines was used before in, e.g., [T8l 1351 fT5l [TOl 
It is known [TB] that, if e -C 5, then MATE has most of the surface area of 
so most pure states are in MATE. This follows from the fact that the average 
of (0|Peq|0) over is equal to ([3]); since (0|Peq|0) cannot exceed 1, it must be 

close to 1 for most 0 (see Section 0] for more detail). It can similarly be shown, see fl?I|) 
below, that most energy eigenstates are in MATE. To be sure, there are states in the 
energy shell which are not in MATE; for example, one could take a tensor product of 
states of two regions having (what look macroscopically like) different temperatures. 

An alternative definition due to Tasaki [15] (in the same direction as mm) , not 
strictly but approximately equivalent and denoted TMATE here, avoids the step of 
rounding off to make the macro observables commute, which may pose substantial dif¬ 
ficulty to carry out in practice. Instead, take Mi,..., Mx to be the macro observables 
before rounding off and coarse graining (mathematically, any self-adjoint operators), let 
Vj = tr(/)™^A0) be the thermal equilibrium value of Mj, and let AM^ be the macro 
resolution of the observable represented by Mj. Using to denote the characteristic 
function of the set A, we define 

Pj = 1[v,-ami,Vi+ami](Mj) (5) 

to be the projection associated with the eigenvalues of Mj that he within the macro 
resolution of the thermal equilibrium value. Then let 

K 

TMATE = njV' e §(^mc) : (01^00) > 1 - . (6) 

j=i 

Note that (0|P00) is the probability of finding, in a quantum measurement of Mj on a 
system in state a value that is AMj-close to Vj. If this probability is > 1 — 5 for at 
least the fraction 1 — rj oi S{J^c) for each j, then TMATE has at least size 1 — Kr] (in 
terms of normalized surface area), which is close to 1 if 77 -C K~^. We note further that 
MATE as in (jT]) can essentially also be written as the right-hand side of ([2]) if the Mj 
are taken again as commuting and coarse-grained on the scale AMj, so that Pj becomes 
the projection to the eigenspace of Mj with the dominant (most degenerate) eigenvalue, 
and J'fTeq is the intersection of these eigenspaces. Indeed, the right-hand side of ([6]) is 
then contained in that of (jTj) with 5 replaced by K5, and that of (jl]) is then contained 
in that of ([H]). 
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3 MITE 


While MATE implies thermal behavior only for macro observables, MITE involves also 
“micro” observables, more precisely, those observables concerning only a region smaller 
than a certain length scale 1. The dehnition of MITE is inspired by canonical typicality^ 
the observation [3 Ea EDI DEI that for any not-too-Iarge subsystem S and most wave 
functions -0 in the energy shell the reduced density matrix of S is close to the 

thermal equilibrium density matrix of S, 


Ab ^ amc 

Ps PS 5 


(7) 


where 

= trsc |0)(0| (8) 

is the reduced density matrix of S obtained by tracing out the complement of S', and 


If S is small enough then 


-me _ amc 

Ps - P ■ 


Ps ~ Ps 


for suitable 0 > 0, where the right-hand side is the partial trace, 

P^P = tTs- , 


(9) 

( 10 ) 

( 11 ) 


of the canonical density matrix 

pW = with Z = ii . (12) 

Zj 

As a consequence, for small S', it does not matter whether one starts from or p^^^ 
(this fact is a version of equivalence of ensembles), and we also have 


p| ^ p)f^. (13) 

We will call flTT]) the canonical or thermal density matrix for S'!! We note that if p^ ~ p™'^ 
for some subsystem S' then the same is true for every smaller subsystem S' contained in 
S' (“subsubsystem property” of ([7])), just by taking another partial trace on both sides 
of the approximate equation p^ p™^. 

The system is said to be in MITE^ (MITE on the length scale S) whenever 0 G 
satisfies ([7]) for every subsystem S corresponding to a spatial region of diameter 
diam(S') < i, i.e.. 


MITE,= f| {0 e §(=^^c) : ~ pr} (14) 

S : diam(5)<£ 

^The density matrix exp{—pHs) with Hs the Hamiltonian of S is sometimes called the canonical 
or thermal density matrix for S'; it agrees with (HU if the interaction between S and its complement 
can be neglected. If the interaction cannot be neglected, then (ED is the correct density matrix to use. 
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with some precise definition of ~ (snch as the trace norm of the difference being smaller 
than a given value). The subsubsystem property implies that every ip in MITE^ lies also 
in MITE^/ for any smaller scale 0 < £' < i. 

MITE is then defined to mean MITE^^ with Iq the largest i small enough to ensure 
that (ITOD holds for every subsystem S with diam(5') < io- As a practical value, for 
example, we may take 

io = 10“^diam(A), (15) 

where A C is the volume of the whole system. Thus a state ip is in MITE if for every 
subsystem of diameter io or smaller, the reduced density matrix is close to the thermal 
equilibrium reduced density matrix. 

Most Ip G S(J^c) lie in MITE. Indeed, canonical typicality (in the sense of pg ~ 
for most Ip) holds for subsystems of up to nearly half the size of A (that is, half the 
degrees of freedom, or square root of the Hilbert space dimension, in practice usually 
half of the volume; see [2911301113] for a discussion). We can choose a moderate number r 
(e.g., r = 8 for cube-shaped A) of overlapping regions S', C A (e.g., also cubes) of nearly 
half the volume so that most ip satisfy ~ p^p for all 1 < f < r simultaneously, and 
so that every region S with diam(S') < | diam(A) is contained in one of the Si. By the 
subsubsystem property, also p^ ~ p™'^ for such regions S, so most ip lie in MITEdiam(A )/4 
and a fortiori in MITE = MITE^^ with io as in (lT5jl . 

A concept along the lines of MITE was used before in, e.g., [32l ESI l36l 125] . 


4 Discussion 

We organize our discussion under a number of subheadings. 

MITE implies MATE .—To discuss this point, it is helpful to first introduce a common 
framework for formulating MITE and MATE. For any observable A, let denote the 

probability distribution defined hj ip ^ over the spectrum of A, 

= (^llB(i)l,*) (16) 

for all sets i? C M (with In (A) the projection to the subspace spanned by the eigenvec¬ 
tors of A with eigenvalue in B). Likewise, let denote the probability distribution 
defined by the micro-canonical ensemble, 

Hy(B) = tr(lB(i)p“) (17) 

for all H C R; p™'^ is the average of p^ with ip taken to be uniformly distributed in 
S(^mc). 

Both MITE and MATE are ultimately of the following form (see also H): For a 
certain family of observables, consider the set of ip ^ for which p^ p™'^ for 

all A G s?/. MATE is obtained by taking .c/ = = {Mi^ ■ ■ ■, Mr}, and MITE^ by 
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taking with the union taken over all regions in A of diameter < I 

and the set of all self-adjoint operators on more precisely, 

£/s = {^0 ® hs'c : Aq self-adjoint on J^s} , (18) 

where I denotes the identity operator and S'^ again the complement of S. Indeed, for 
A = Aq<^ isc and G MITE^, 

,,y(S) = tr(lB(i)p-) = tr(lB(i„)pr) « tr(lB(i„)pp = (19) 

for all -B C M by (IT^ . 

From this perspective it is obvious that MITE implies MATE when based on rea¬ 
sonable choices: Suppose that 

L<io, (20) 

where L is the length scale of the macro observables—the diameter of the cells Aj on 
which the macro observables Mj were dehned at the beginning of Section [2l that is, 
suppose that ([7]) holds at least up to the length scale of the macro observables. This is 
commonly the case; e.g., for a cubic meter of gas at room conditions, we can realistically 
take L ^ 10“'^ m and io ~ 10“^ m. As a consequence of fl20l) . jz^iate C ^mite, so if 
-0 G MITE then ^jJ G MATE. 

ETH. —To come back to the eigenstate thermalization hypothesis (ETH), it comes in 
two variants: MATE-ETH and a more rehned version MITE-ETH, according to whether 
the energy eigenstates are required to be in MATE or MITE. It is MITE-ETH that fails 
dramatically in some MBL systems, according to the Endings of El ESI DS]; there it 
is shown for certain MBL systems that a substantial fraction of the energy eigenstates 
(in a micro-canonical energy interval), or even all of them, he outside of MITE. At 
the same time, it is easy to see that for every macroscopic quantum system (MBL or 
not), MATE-ETH must be almost satisfied, in the sense that most energy eigenstates 
\n) G S(J^c) are in MATE: Assuming that ([3]) holds with £ -C h, we obtain, writing 
D = dimJ^jnc, that 

1 ^ ^ 1 ^ 

— ^(n|Aq|n) = — tr(Peq) = 1(21) 

n=l 

and since {n\Peq\n) cannot exceed 1, most of these terms must be close to 1. 

If MATE-ETH holds strictly, i.e., if all energy eigenstates in J^rac are in MATE, then 
every state -0 G S(JifAc) will sooner or later reach MATE and spend most of the time 
in MATE in the long run. That is because [15], writing /(f) = lini'r^oo ^ Jq f{t) dt for 
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time averages, \n) for the energy eigenstate with eigenvalue Em and = e 


= ^(^|n) e*®"‘(n|Peq|?^')e {n'\'il^) (22) 

n^n' 

= '^\{P\n)\‘^{n\Peq\n) (23) 

n 

> j;|(^,|n)|'^(l-i) (24) 

n 

= 1-5, (25) 


provided H is non-degenerate, i.e., ^ E^' for n ^ n' (using e*®* = 1 if P = 0 and 

= 0 otherwise)!^ Since its time average is close to 1, ('0t|Feq|'0t) must be close to 1 for 
most t in the long run. 

It follows from this that systems with MBL for which the transport coefficients van¬ 
ish, so that an initial state ip with a non-uniform temperature will remain so indehnitely, 
cannot have all of its energy eigenfunctions in MATE. Since most energy eigenstates are 
in MATE, such ip must be a superposition of predominantly those rare eigenstates that 
are not in MATE. 

This leads to the question whether there are macroscopic systems for which all 
energy eigenstates are in MATE—i.e., whether MATE-ETH ever strictly holds. It is 
known that this is so for a random Hamiltonian whose eigenbasis is uniformly chosen 
among all orthonormal bases [15] ; see also uni HI]. Some numerical evidence [20] points 
to the existence of systems with realistic interactions for which all energy eigenstates 
are in MITE and thus also in MATE. 

For MITE-ETH, there are several results jSS] HHl E3] , all of which assume that the 
Hamiltonian is non-degenerate and has non-degenerate energy gaps, i.e.. 


Em Eyi ^ Em' E^i' unless 


either m = m! and n = n' 
or m = n and m' = n', 


(26) 


a condition that is generically fulhlled. We note here two results, the hrst of which [23] 
asserts that if all energy eigenstates in are in MITE, then most ip G will 

sooner or later reach MITE and spend most of the time in MITE in the long run. More 
precisely, those ip will behave this way for which the effective number of signihcantly 
participating energy eigenstates is much larger than dim for any small S. The second 
result [35] shows that all (rather than most) ip will ultimately reach MITE and stay 
there most of the time, under two assumptions, hrst again that all energy eigenstates are 
in MITE, and second Srednicki ’s [nmg extension of the ETH to off-diagonal elements, 
i.e., that for A G . 2 / (here, ssf = Us-^s as in flTSll h 

(m|A|n) 0 for m 7 ^ n (27) 

^In fact, the assumption of non-degeneracy can be dropped: If we number the eigenvalues as 
with En En' for n ^ n' and let \n) denote the normalized projection of ip to the eigenspace of E^, 
then the calculation (|22p -(l25 1) still applies. 






(see also H). Indeed, a calculation using fl26|) shows that 


\{m\A\n)\‘^ \{n\'ip)\‘^, (28) 

^ / m^n 

and if |(m|A|n)| < e -C 1 for all m ^ n, then the time variance (l28il is smaller than e^. 
If all \n) are in MITE, a calculation similar to (l2^ - (]25ll shows that 

(V>t|i|V't)~tr(/)“^i). (29) 

It follows that, for most t in the long run, ~ A) for all A E Us^s (in 

particular for projections), so tfjt E MITE for most t in the long run. 

Mixed states. —Once we have the notions of MITE and MATE for pure states 'll), 
they are easily generalized to mixed states p\ MATE occurs if tr(.Peq p) > 1 — 5, and 
MITE if ps ~ for all subsystems S dehned by spatial regions of diameter < Iq. Note 
that neither MATE nor MITE requires that p be close to or 

Thermal equilibrium in classical mechanics. —Only one of the two notions MITE and 
MATE can be satished for pure states in classical mechanics, namely MATE. That is 
because a “pure state” corresponds in classical mechanics to a point X in phase space, 
while a “mixed state” corresponds to a probability distribution over phase space. Since 
X specihes the positions and momenta of all particles, it also provides a pure state for 
any subsystem. In contrast, in quantum mechanics pg can be mixed, and in fact is mixed 
except for product states. So in classical mechanics it is never true for a system in a pure 
state that a subsystem S could have a state close to a thermodynamic ensemble such 
as the marginal (obtained by integrating out the variables not belonging to S) of the 
micro-canonical distribution (i.e., uniform over the energy shell) or the canonical one for 
the whole system. In contrast, MATE is analogous to Boltzmann’s [H [9l [Til [22] notion 
of thermal equilibrium for a closed classical system, based on a partition of phase space 
into macro states T^. (Note that there is no difference between MATE and TMATE 
classically, as all observables commute.) 

Abstract MITE. —A natural mathematical generalization that is often interesting to 
consider is based on dropping the idea that S corresponds to a region in 3-space and 
regarding S as an abstract subsystem dehned by any splitting of Hilbert space into a 
tensor product, 

^mc C ^5 0 ^5. , (30) 

where S''’ can be thought of as just an index for another Hilbert space. For example, S 
may comprise the spin degrees of freedom and S'"’ the position degrees of freedom, or S 
may comprise the oxygen atoms and S''’ all other atoms in the system. Given a list of 
subsystems Si,..., S'^, one can dehne 

r 

MITEs, .s, = 6 S(JC,e) : pt. « «“} . (31) 

i=l 
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Canonical typicality implies that, if r is not too large and each Si is not too large (a snf- 
hcient condition is dim.;?^, -C \/dim J^c), then most 'ijj E S(J^c) are in 
see Theorem 1 in [2S1 ED] and Proposition 1 in [TB] for a precise and qnantitative state¬ 
ment of this fact. 

One can also consider the set MITEmost comprising those G §(J^c) for which 
pt ~ holds for most abstract snbsystems S of dimension < cIq. If do ^ \/dim J^rac, 
then also MITEmost has most of the snrface area of S(J^c) [IS]- On the other hand, 
given any pnre state 'ijj E ~ cannot hold for all abstract snbsystems S 

of dimension < do simnltaneonsly [25 |T3] . 

5 Conclusions 

Perhaps the most snrprising aspect of the sitnation is that the varions criteria for thermal 
eqnilibrinm of pnre states proposed in the literatnre fall into two gronps that differ 
snbstantially in how mnch they demand. 

Argnably, the essence of thermal eqnilibrinm is what characterizes it in thermody¬ 
namics: that a system appears stationary on the macro level, and that temperatnre 
and all chemical potentials are spatially nniform. This corresponds to MATE, which 
may therefore be regarded as the direct expression of thermal eqnilibrinm. On the other 
hand, since MITE is the stronger statement, and since it is usnally trne that macroscopic 
qnantnm systems approach MITE (MBL systems being an exception), it is natnral to 
consider MITE, and it wonld seem artificial to not regard it as a new kind of thermal 
eqnilibrinm property emerging from qnantnm entanglement. 

Acknowledgments. J. L. Lebowitz was snpported in part by the National Science Fonn- 
dation [grant DMR1104500]. 
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